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LINEARIZED THEORY OF STAGNATION REGION SKIN FRICTION 


AND MASS TRANSFER AT HYPERSONIC SPEEDS 
Kenneth K. Yoshikawa 
Ames Research Center 


SUMMARY 


A linearized theory for stagnation region heat and mass transfer has been 
extended to include the shear stress. The method considers the momentum and 
energy equations which are uncoupled first and then linearized by the intro- 
duction of the effective values for the thermal and transport properties. 
Solutions of these equations are applicable at very high speeds, where conduc- 
tion and shear stress in the shock layer and the interaction between the 
heated gas behind the shock wave and the injected or ablated gases must be 
considered. An explicit closed-form first order solution has been obtained 
for the compressible, viscous shock-layer momentum equation with mass 
addition or suction. 

Results of this analysis show that the effect of gas injection on shear 
gradient (skin friction) near the stagnation point strongly depends on the 
Stanton number without mass addition, wall temperature, relative injectant 
molecular weights, and the ratio of the averaged Peclet number between the 
free-stream and injected gases. The solution is substantiated by comparisons 
with exact solutions. 


INTRODUCTION 


Significant reductions in heat transfer and skin friction in the presence 
of foreign gas injection have been reported in a number of papers (refs. 1-7). 
These reductions are extremely important in the design of heat and shear 
resistant vehicles for very high speeds. Theoretical approaches to analyze 
the effect of mass addition on heat transfer and shear stress employ either 
(1) boundary- layer theory, or (2) shock-layer theory. The results may differ 
significantly, depending upon the flight conditions. For very high speed 
flight, shock- layer theory is essential for analyzing properly the problems of 
heat transfer and shear stress in the presence of mass addition. 

Although the skin friction vanishes at the stagnation point of a blunt 
body, the ratio of shearing stress with mass addition to that without mass 
addition has a nonzero value. This ratio is important because it is related 
to the reduction in skin friction along the body surface, that is, the larger 
the reduction in gradient at the stagnation point, the smaller the 
skin- friction along the body surface. 



The primary purpose o£ this report is, therefore, to find a simple 
solution, for the reduced shear, by extending the method of reference 1 to the 
shock-layer momentum equation. To find a solution the momentiom equation is 
decoupled from the other conservation equations (continuity and energy equa- 
tions) and linearized by means of an effective Reynolds number. A closed-form 
solution demonstrates the dependence of shear stress on the appropriate 
physical parameters. The method is substantiated by comparison with a number 
of numerical solutions of the exact equations. 


ANALYSIS 


A simplified momentum equation for the viscous, compressible flow in the 
neighborhood of a stagnation point has been considered. The reduction of the 
shear stress at the wall by transpiration or ablation gases emanating in 
hypervelocity flight will be emphasized. The present analysis considers the 
momentum and energy equations that will be decoupled by approximating the mass 
flow term in the equations . The equations are then linearized by the intro- 
duction of the effective values for the thermal and transport properties of 
the gases. To provide a physical insight into the parameters, a first-order 
solution is derived from the linearized equations for the case of no heat 
conduction and shear stress just behind the shock wave. 

Finally, the primary effect of heat conduction and shear stress 
(postcooling effect) behind the shock wave on the flow parameters is 
determined by solving the linearized equations. 

Although the shear stress at the stagnation point is zero, the ratio of 
shear stress with gas injection to that without gas injection is finite. In 
what follows, this ratio is designated by $ and represents the reduction in 
shear gradient at the stagnation point. It may also be considered as a mea- 
sure of the reduction of shear stress (skin friction) in the vicinity of the 
s t agnat i on po in t . 


Geometry and Assiomptions 


Figure 1 shows the geometry used in the analysis, 
as follows: 


The assumptions are 


1. Steady axisymmetric stagnation-point 
flow 

2. Thermodynamic equilibrium 

3. No chemical reaction between free- 
stream and injection gases 

4. Laminar continuum flow 

5. Thin shock layer, L/R << 1 

6. Small blowing rates, Pw^w/Pcx>^oo ^ 

7. Thermal diffusion and diffusion- 
thermal neglected 
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Basic Equations and Boundary Conditions 
The basic equations used for the present analysis are: 
Continuity equation: 


(pur) + (Qpvr) = 0 


( 1 ) 


where 


Q= 1 -I 


Momentum equations : 


3u ^ 3u 


9y 


0 



( 2 ) 


Energy equation: 


Since the major portion of this study is concerned with the reduction of 
shearing stress in the vicinity of the stagnation point, the momentum and 
enery equations for the shock-layer flow are further simplified, for 
9p/8x - -p^Ug (3u/8x) g , pu(3j/3x) << pv(3j /9y), and (1/2) (u^ + v^) << h; 


-Qpv 


3u ^3 



pu 


3u 

3x 


- p u. 

s 6 



(4a) 



(4b) 


The boundary conditions immediately behind the shock wave are altered by heat 
conduction and shear stress in the shock layer. The momentum and energy 
balances behind the shock wave may be written as 
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,, X 1 / 9u\ 

^ ^ “ p J. V ^y)s 

^ P«.Voo \Cp 9y/s 


where H is the free-stream total enthalpy. 

Equations (4) are solved independently by the introduction o£ an equation 
for the normalized mass-flow distribution which is relatively independent of 
heat transfer. The mass-flow equation taken from reference 1 is 

X = = Xw - Xw’n - (1 - x„' + Xw)n^ (6) 

^OO 00 

where n = y/L is the distance from the wall normalized by the shock-standoff 
distance, and Xw^ is the absolute value of 3x/9n near the stagnation point 
as shown in reference 1. 


Transformed Equations 

From equations (4) and (6) , the dimensionless momentum and energy 
equations for the shock layer after some minor manipulation and approximation 
become (see appendix A) 


d / 1 du \ 
dn \ Re dn / 


X ^ ' '‘a ^ PCa) 


A. (L . y Al - 

dn \3 dn/ ^ dn 


where the dimensionless quantities are defined as 


— _ u u 

■ u " v::s7r 


H - 


[CusX'/Xsn^- 2e(l-e)] 


F(n) = 


PooV^^L 


p V L 

OO 00 


Q = 1 + kae n 


1 + /(8/3)e - e 


The physical quantities used in equations C7) are described in detail in 
appendixes A and C. 
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The boundary conditions taking the blowing from the wall into account and 
using the modified Rankine-Hugoniot condition across the shock wave are, from 
equations (5) : 

at n = 0 


at n = 1 


X = x„ 


u = 0 


h 




C8a) 


(8b) 


Equations (7) are nonlinear second-order (ordinary) differential equations 
with nonlinear boundary conditions (eqs. (8)) and can only be solved 
numerically. 


Linearized Solutions 

Linearized solutions to the momentum and energy equations can be obtained 
explicitly when the Peclet number 8 and Reynolds number Re are replaced by 
values considered constant for the specified boundary temperature (see ref. 1 
and appendix B) ; for example, for a cold wall (T^^ = 0), 


8 = 6 * 
Re = Re* 


( 9 ) 


The linearized solutions become 


u 


h - 


‘w 




ClOa) 


£ ^ceVRe*)G(m) 


(10b) 
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where 


c = 1 



GCni) 


Re 


■i™ ■&). £ 


(lOc) 



dn^ 



G(n) = 


J_eG(l) 

Re* 



Re* f X drij 

Jo 




Note that Uq, the complementary solution to equation (7a) , has the same form 
as equation (10b) with the constant g* replaced by Re*. For the general 
wall condition (T^ 0) the quantities marked by an asterisk in equations (10) 

will simply be replaced by double asterisks. 

Since detailed analysis of the linearized solution to the energy equation 
and the heat reduction function ip have been described in reference 1, only 
solutions to the momentum equation and the reduction in the shear gradient 
will be considered in the following sections. 


First-Order Solution Without Postcooling 

It is shown in reference 1 that for a large Peclet number (g >> 1) the 
boundary condition for the flow enthalpy given in equation ^8b) is not 
affected by conduction behind the shock wave and hg = 1 + hj^. Similarly, the 
boundary condition at the shock for the velocity in equation (8b) will not be 
significantly affected by conduction or shear stress for high Reynolds numbers 
(Re >> 1), and Us == 1. It follows that one can derive a simple first-order 
solution for equation (10a) by taking only the leading term from the velocity 
distribution 


u - ciUj. + (1 - ci)n + . . . (11a) 

which satisfies boundary conditions .il(0) = 0 and u(l) = 1 (ci == c) . Since 
the velocity gradient near the shock wave is relatively unchanged by mass 
addition (i.e., c is insensitive to blowing as shown by numerical solutions 
of eq. 10(d). See also ref. 3), the constant ci can be evaluated simply by 
integrating equation (7a) from n = 0 to n = 1 for = 0. Thus, 


6 


^ I u” dn, - I XU' dn = ka I FCn)dn 

Jo Jo Jo 


1 


(11b) 


Substituting equation (11a) into equation (11b) and using the relation 

•1 

- XUc')dn = 0 


f(. 


one can calculate the constant ci as 


Cl ^ 1 - 


2 + Xw' 


- {^) * (^) - 


for e << 1 (Q = 1) (11c) 


Note that is independent of Re*. 

For Re* >> 1 and X^ << 1, the complementary solution (eq. (10c)), will 
become (from ref. 1) 

Ug erf /-j (ReXw')* n - a + Ug erf a 

uc = =! (lid) 

1 + erf a 


The velocity gradient from equations (11a) and (lid) becomes 



Cl 



Uc 



1 + erf a 


+ (1 


Cl) 


( 12 ) 


where 




a 



(Re/x„') 



Xw 


Equations for shear reduction- The reduction in shear gradient with mass 
addition can be obtained in terms of the blowing parameter B by dividing 
equation (12) by the same equation for B = 0. Thus, from equation (11a), the 
shear reduction can be written 



(13a) 
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where the first term in the numerator and the denominator, which is a part of 
the con^lementary solution, can be derived by a method exactly analogous to 
the one used for the heat-transfer equation described in reference 1; that is. 



(13b) 


The blowing parameter and effective Reynolds number (see appendix B) are 
given by 


where 


B = Xw/(St)o 
Re* = B*/(Pr)„ 

(St)o = /(2 /tt) (Xw’/eo)* 


(14) 


Dividing equation (13a) by the first term in the numerator, and incorporating 
equations (12), (13b), (14) and appendixes C and D, one will obtain the 
following result: 


(<|>rx,-r=L= 

, - I- 


(^^/cp)w 1 

(k/Co)* L/Lo 


, (k/cp)w 

^ ^ Tk^^ 


(15a) 


where 


A. 






W- 


1 L 

” T \ 


TT Lq (Pr) 


w 


IL >B/Bt 


1 + erf 


1 L 


/TX 

/ ir Lq 


(P^)w 


IT Lo (Pr) 


w 


1 + erf 


/TX 

V ^ Lo 


Lo (Pr) 


w 


for B 


> (15b) 


for B > Bf 
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♦ 


* 


In equations (15) (Pr)^ is the Prandtl number evaluated at the wall condition. 
Note that the function (f) (eq. (15b)) is ide ntical to the heat-reduction func- 
tion ijj, if the blowing parameter B// (Pr)^ is replaced by B. As shown in 
reference 1, the values of Wf given in. equation (15c) are determined from 
either theoretical calculations or experimental data for the heat-transfer 
parameter by 



(15d) 


For a first-order approximation, the factor Nf 
as the main stream) 


is found to be (for air 


Nf ^ 

5 

9 

for 

Nf 

1 

for 


monatomic gas 

diatomic and some polyatomic gases^ 


Derivations of these values have been discussed in detail in reference 1. 

^N£ - 5/9 for Freon gases (polymolecular), (data deduced from ref. 5). 
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shear reduction for identical-gas infection- As a special case, 
equations (15) can be further simplified to 


$ 


+ Y 


(k/cp)w 1 -| 

(k/Cp)* L/Lo 


(k/cp)„ 

L ^ ^ ^ (k/Cp)* 


(16a) 


where 




1 L 1 2 

g TT Lo (Pr)w 



1 + erf 


1 Jl. 

TT Lo 


1 

CPr)„ 


B 


for all B (16b) 


It will be shown later that equations (15) and (16) are in satisfactory 
agreement with the numerical solutions of references 3 and 4. 

When L/Lq = 1, = 0 and y ~ 1 j equations (16a) may be expanded in 

a series to obtain 


Cl^/cp)w /(Pr)„ 

(k/ Cp) 


(16c) 


Effect of Heat Conduction and Shear Stress 
Behind the Shock Wave 

The enthalpy and tangential velocity gradient immediately behind a shock 
wave (assuming the strong shock relation) will be altered significantly when 
heat conduction and shear stress become major factors in the energy and momen- 
tum balances within the shock layer. As in reference 1, this phenomenon is 
designated here as the "postcooling effect." When this phenomenon occurs, the 
first-order solution does not apply, and numerical solutions of the linearized 
equations (eqs . (10)) incorporating the effective Peclet number (or Reynolds 
number) from table 1^ must be obtained. Preliminary effects of postcooling 
on the shear gradient reduction in the vicinity of the stagnation point will 
be discussed in the section "Results and Discussion." 
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^See appendix B. 



RESULTS AND DISCUSSION 


Calculations o£ the shear reduction at the stagnation point including 
foreign gas injection and suction are presented in this section for atmo- 
spheric flight conditions where the first-order solution applies. The post- 
cooling effect on the shear stress at relatively low Reynolds numbers is also 
examined. 


Results Without Postcooling 


Shear veduot-ion - ■ident'Loat gas injections- A comparison of the present 
first-order solution (eqs . (16a) and (16b) for p^ = 1 atm, Tg = 15,000° K, 
and R = 30 cm (1 ft) with the exact numerical solutions of Howe and Sheaffer^ 
(ref. 3) is presented in figure 2 for identical gas injection. The present 



B 


Figure 2.— Effect of gas injection on shear 
reduction — air injection (pg = 1 atm, 

Ts ^ 1 5,000^^ K, and ^ 1 ,000° K). 



Figure 3.— Effect of Stanton number on shear 
reduction — air injection (pg = 1 atm, 

Ts ^ 1 5,000° K, and ^ 1 ,000° K). 


results show satisfactory agreement with 
the numerical solutions. Also presented 
is the solution obtained with the linear 
equation (16c) , which is very close to 
the first-order solution for B < 1. 

Effect of Stanton ntunber- The 
present solutions show that the shear 
reduction function as for heat 
reduction (ref. 1), depends on the 
stagnation-point Stanton number (St)^, 
which is a function of the product of 
free-stream density and body radius 
(ref. 1) . As the stagnation-point 
Stanton number (St)o increases (i.e., 
the product of free-stream density and 
body radius decreases) , the mass trans- 
fer becomes more effective in reducing 
the shear stress, as illustrated in 
Fig. 3. The present solutions and those 
of Howe and Sheaf fer agree closely. A 
single correlation curve for the shock- 
layer theory can be obtained only if the 
Stanton number (St)^ remains constant 
and T^/Ts - 0. The effect of the 
Stanton number (or p^R) on shear stress 
(as well as heat transfer) is due essen- 
tially to the changes in the shock- layer 
thickness, which are accompanied by 
changes in the velocity gradient at the 
stagnation point. Experimental data on 
heat transfer and pressure changes due 
to the mass injection support the 
present conclusion (ref. 5) . 


^By courtesy of Howe and Sheaffer since numerical results for shear 
reduction were not presented in reference 3. 
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The effect of transport properties on shear reduction is also recognized; 
that is, (in fig. 3) the shear reduction curve for p^ = 10 atm and R = 3 cm 
(0.1 ft) is separated from that for pg = 1 atm and R = 30 cm (1 ft) even 
though the Stanton number and wall temperature for both cases are about the 
same. This is due to the pressure dependence of the thermal and transport 
property parameter (k/Cp),^/ (k/cp) * . The heat-reduction function however, 
has a single correlation curve when the Stanton number is kept constant as 
demonstrated in reference 1. 


PRESENT METHOD 


Effect of wait tempevaticre- In contrast to the heat reduction, the shear 
reduction is sensitive to the wall temperature. This dependency can be recog- 
nized immediately from the result of the first-order solution since the quan- 
tity (k/cp)^ in equation (15a) depends on wall temperature. The effect of 
wall temperature on shear reduction is shown in figure 4 for air injection. 

Generally, the higher the wall tempera- 
ture, the less the shear reduction. The 
present methods are compared with the 
numerical calculation of reference 4, and 
the results show satisfactory agreement 
for a wide range of wall temperatures. 

To make this comparison, the parameter 
(St)Q, in equation (C8) , is main- 
tained at 0.15, the value used for the 
heat transfer case. 


Pg, atm 
. 1-10 
10 
10 


ref 4 
\ T^.”K 

0-3 300 

3/4 3,000 

I 300 



Figure 4 Effect of wall temperature on 
shear reduction — air injection 
(VT/F (St)o=- 0.015). 


■ PRESENT METHOD 
(eqs(i5)) 


Shear reduction - foreign gas 
injection- Injecting gases different from 
the free-stream gas affects the shear 
reduction in the manner shown in fig- 
ure 5. The present results, derived by 

means of equations (_15) , for Ps = 1 atm, 

Ts = 15,000° K, /ITT (St)o = 0.15, and 
for Tw/Tg = 0, are compared with the 
boundary-layer solutions from reference 4. 
Satisfactory agreement is obtained over 
the entire blowing range including suc- 
tion (B < 0) . As can be recognized from 
equations (15) , shear reduction for a 
cold wall (Tw/Ts = 0) depends primarily 
on Prandtl number at the wall and on the 
blowing parameter B. The injection of 
light gases is more effective in reducing 
the rate of the shear stress near the 

Figures.— Effect of foreign gas injection on shear stagnation point. Note that no combus- 
reduction (pj = 1 atm, Ts = 15,000° K, and tion is assumed to take place when 
T\v/Ts — 0). hydrogen is injected. 
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Effect of Postcooling on Shear Reduction 


The effect of postcooling on the stagnation-point flow is (1) to reduce 
the driving enthalpy which changes the transport properties, and (2) to reduce 
the effective tangential velocity due to the shear stress in the shock layer, 
which lowers the shear stress at the body surface. As shown in reference 1, 
the enthalpy behind the shock wave varies with Stanton number and blowing rate. 
The present analysis gives a similar result for tangential velocity as shown 

in figure 6. The results were obtained 
from numerical solutions of equation (10) 
since the first-order solution neglects 
the postcooling effects.^ The immediate 
effect of the shear stress on the velocity 
ratio behind the shock wave is noted. 

This is due to the nonzero value of the 
velocity gradient (eq. (11a)) and thus to 
the finite value of shearing stress ratio 
at the shock wave. Since the velocity 
(also enthalpy) varies with Stanton number 
and blowing rate, the shear reduction for 
air injection is altered in the manner 
shown in figure 7. For a small Stanton 
number the shear reduction is the same as 
obtained from the first-order solution 
(e.g., fig. 3). The shear reduction 
becomes less effective, however, as the 
Stanton number increases beyond the crit- 
ical number, (St)^^ - 0,32 (see ref. 1), 
and, finally, $ approaches the free 
molecular limit of unity as (St)^ 1. 

Similar results are shown in reference 1 
for the heat reduction. 



Figure 6 — Effect of Stanton number 
on tangential velocity behind shock 
wave — air injection (pg = 1 atm, 

Ts ^ 15,000'" K, and Tw/Tg 0). 



The critical Stanton number, at which 
the effect of postcooling becomes impor- 
tant, can be defined in terms of the phys- 
ical parameters through the effective 
Peclet number 3* (see eq. (7c)). From 
figure 10(a) of reference 1, the effective 
Peclet number is found to be approximately 
3* - 5 for (St) 0 ]^ “ ^ .32. Thus, one can 
show, in an actual flight case where 

- 15 km/sec, that the effect of postcooling on convective heat and shear 
stress becomes important when the product of the free-stream density and body 
radius becomes p^R < 3xl0""^ cm, which corresponds to a flight altitude of 
58 km for a body radius of R < 1 cm, or a flight altitude of 84 km for a 
body radius of R = 30 cm (1 ft) . 


Figure 7.— Effect of Stanton number on shear 
reduction — air injection (pg = 1 atm, 

Ts ^ 15,000° K, and Tw/Ts ^ 0). 


^The computing program has been prepared by Mrs. Sarah Rogallo, NASA, 
Ames Research Center. 
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CONCLUDING REMARKS 


Closed- form solutions have been obtained for the reduction of the shear 
gradient normal to the surface in the vicinity of the stagnation point of a 
body with mass transfer. The shock-layer momentum equation was decoupled from 
the continuity and energy equations by approximating the mass- flow distribu- 
tion along the stagnation stream line. The momentum equation is then linear- 
ized by the introduction of an effective Reynolds number which is a constant 
for a given boundary temperature. The linearized solution indicates conclu- 
sions similar to those for the heat transfer for the reduction of the shear 
gradient at the stagnation point by mass addition, which represents the 
direct reduction of skin friction in the vicinity of the stagnation point. 

1. The degree of reduction depends on the mass-addition rate and the 
Stanton number for no mass addition. For a given mass-addition rate, the 
shear- reduct ion rate increases with increasing Stanton number until the 
Stanton number reaches the critical value of about 0.32. The shear reduction 
becomes less effective as the Stanton number increases beyond the critical 
number and approaches the free molecule flow limit of unity. 

2. The shear reduction depends on wall temperature; the higher the wall 
temperature, the less the reduction in shear stress. 

3. The injection of gases of lower molecular weight is most effective 
for the shear reduction. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, June 11, 1970 


APPENDIX A 


TANGENTIAL MASS FLOW GRADIENT 


The relation between the tangential and normal mass flow gradient in the 
vicinity of the stagnation point can be derived as follows: 

The continuity equation (eq. (1)), with the effect of shock-layer 
thickness taken into account, is 

(pur) + ^ CQpvr) = 0 

where Q = 1 + y/R and r = Qx. If it is assumed that pu = [9(pu)/3x]x and 
pv [9(pv)/9y]y, the tangential mass-flow gradient and momentum become 

^ - (I Q - l) ^ tAI 

- (I7) “ ^ 

To allow for the velocity gradient at the shock wave, equation (Al) is 
modified as 


9(pu) 

9x 



(A3) 


where the factor c^ is determined from the boundary condition for the tan- 
gential mass- flow gradient at the shock wave. 


9(pu) 

3x 



(A4) 


Thus , from equations (A3) and (A4) , c^;* is given by 


Cr " 




(AS) 


For hypersonic flow where e << 1 and << 1 it follows, from equation (A3), 
that 


15 




Thus , equation (A2) incorporating equation (A6) can now be written as 



which is the result introduced in equation (7a) . 


(A7) 
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APPENDIX B 


THE EFFECTIVE THERMAL AND TRANSPORT PROPERTY PARAMETERS 


Values for the effective Peclet number and Reynolds number which will be 
used in the linearized solution (ec[s . (10) and (15)) are derived in this 
appendix. 


THE EFFECTIVE PECLET NUMBER 


Let 6 1 be defined as 


where q = [ (k/Cp) (dh/dy) ] and also 6^ to be 

(dh/dy) dy 

At = •'O 

- (dh/dy) „ 


(Bl) 


(B2) 


It follows, from equations (Bl) and (B2) , that 

■ (5)'(a (i) 

where 

A y , /o" ^ ^ fo"" 

V^p/ (dh/dy) dy dh 






for h„ = 0 


(B3) 


Note that the averaged value (k/Cp)* requires no knowledge of the solution. 
Now, assume that the variable (k/Cp) can be replaced by some appropriate con- 
stant value, say (k/cp) = X; it then follows, from equations (Bl) , (B2) , and 
(B3) , that 


and 



\ 





(B4) 
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■ Ml 




CB6) 


Table 1 presents some numerical calculations of the effective thermal and 
transport property parameter which are used in the first-order and general 
solutions. Table 1 seems to provide better results and has a more logical 
basis than that given in reference 1. 


TABLE 1.- THERMAL TRANSPORT PROPERTIES, (k/cp)*/(k/cp\ef 



0 0.0000000 0.0000000 0.0000000 

1000 .6126295 .6126295 .6126295 

2000 .9499233 .9502330 .9520912 

3000 1.2197483 1.2552291 1.3579450 

4000 1.5759500 1.7378318 1.8327561 

5000 1.9541369 1.9647090 1.9875224 

6000 2.1406367 2.2481930 2.4740788 

7000 2.4545306 2.7457312 2.9321786 

8000 2.9100277 3.1459877 3.0117004 

9000 3.3108144 3.2662349 3.0849595 

10000 3.5217905 3.3439494 3.4486317 

11000 3.6278065 3.6939754 4.1170344 

12000 3.7265962 4.2138305 5.0633505 

13000 3.9544243 4.8933558 6.1638561 

14000 4.4721040 5.7415978 7.1545277 

15000 5.1004932 6.6690641 7.8287131 


10'^ 

10’^ 

10"^ 

10-^ 

0.0000000 

0.0000000 

0.0000000 

0.0000000 

.6126295 

.6126295 

.6126295 

.6126295 

.9559015 

.9728052 

1.0208307 

1.1343960 

1.4889271 

1.5654652 

1.8207278 

1.5203098 

1.7297006 

1.7252176 

1.9953843 

1.9555636 

2.0680217 

2.3124995 

2.5002307 

2.4097264 

2.6755975 

2.6202073 

2,4553741 

2.4686934 

2.7763997 

2.6461656 

2.7342026 

3.0686967 

2.8403820 

3.0464870 

3.6058875 

4.6225201 

3.3162889 

3.9997432 

5.1343087 

5.6699780 

4.2053297 

5.4051008 

6.0928483 

5.7321583 

5.4612288 

6.4812224 

6.2490610 

5.6746416 

6.6742440 

6.8617359 

6.2177217 

5.6041260 

7.4202914 

6.9172425 

6.1587323 

5.5313316 

7.7161356 

6.8953037 

6.0941142 

5.4617239 

7.8154403 

6.8470945 

6.0360700 

5.4029145 


\/cp =ju/Pr ; (k/cp)jgf = 36.9X10 * Ibm/sec-ft or = 550X10 ^ gm/sec-cm; T^gf “ 1.000° K. 
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THE EFFECTIVE REYNOLDS NUMBER 


Accurate calculations of heat transfer to the surface can be obtained if 
the average Peclet number 3** is used for general wall temperatures in the 
energy equation. Similarly, the effective Reynolds number Re** may be 
introduced into the momentum equation where Re** is related to 3** by 


(Pr)** 

where (Pr)** is approximated by 


for ~ ^ 0 
As 


(Pr)** = (Pr) 


(B7) 


For T^/Tg - l,(Pr)** becomes the Prandtl number for the wall temperature, and 
n = 1. For Tvv^/Tg - 0 the results of reference 2 (ch. 4) yields 
Re* - 3*/(Pr)^^'^, and it can be shown that (Pr)^^^ ~ C^^)w since the Prandtl 
number is near unity for the lower wall temperature. Thus, the effective 
Reynolds number covering the intermediate wall temperature case is 
approximated simply by 


(Re)** 


g * * 


(B8) 
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APPENDIX C 


EQUATIONS FOR THE PHYSICAL PARAMETERS 
BOUNDARY- LAYER THICKNESS 


An important physical dimension in the viscous layer is the 
characteristic thickness 6+, arbitrarily defined in equation (B2) as 

^ (dh/dn)dn 

(dh/dn)^ 

Since from reference 1, 




it follows that 



(Cl) 


The thermal boundary- layer thickness for 


0 is from reference 1, 



(C2) 


Thiis , 


(^l “ / (6oV)* 


CXw')* 


(C3) 


Similarly, the characteristic thickness corresponding to the momentum equation 
(from the homogeneous solution, eq. (12)) may be written as 



(C4) 
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The velocity boundary- layer thickness is, thus 




■ (r) = / 


2tt 


(ReoXw ’ ) 


(C5) 


BASIC FLOW PARAMETERS 


The approximate relations for the basic flow parameters taken from 
reference 1 are : 

Mass- flew gradient for 1^/1^ - 0: 

(Xw')* ^ ka /8e(l - e) (C6) 

(see eq. (7c) for k^) . 

Mass- flow gradient for T^/T^ 0 : 

(Xw')** " k;,(Xw')* (C7) 


where the factor = (k/Cp) */ (k/cp) ** depends on wall ten 5 )erature (ref. 1) . 

Shook-tayev thickness : 

^ 1 ^ (mI) kw 

or, from equation (13) , 

1 (St)oB (C8) 

where k^, depends on the body shape (k^^ = 1.0 for an axisymmetric body). 

The effective htaioing ’parametev for the hoimdaxy-layeT btowoff : 

Blowoff is the situation where the increase in the shock-layer thickness 
is equal to the original boundary- layer thickness. From equations (C3) and 
(C8), 


TT/E(Mf/M^) 

cx„')* 


(C9) 
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Similarly, from equation (C5) 


B 


t 


Tr/(Pr)„eCM£7M«.) 

Cxw')* 


/(Pr)w B* 


Equation (CIO) is given in equation (15c) . 


(CIO) 
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APPENDIX D 

FIRST-ORDER SOLUTION TO THE MOMENTUM EQUATION 


A closed form solution to the linearized equation (7a) is given by 
equation (10a). An explicit closed form solution, is obtained by taking only 
the leading term from the linearized solution. The shear reduction function 
$ can be written, from equations (13), as 


$ 


T (uc'/Re)w,B 1 ^ 

n - ci>| 

M * V 1 

j_(uc'/Re)^^B»o_ 

\ Cl J 

'^®w,B V«c'/Re4^B=0 1 


1 

( 1 y ( 

^ \ Cl ) 

R®w,B=o 

\uc’/Re/„^B=o } 


(Dl) 


The first term in the numerator of equation (Dl) takes exactly the same form 
as the heat red uction function except that the blowing parameter B is 
replaced by B//(Pr)Y^, 


(Uc'/Re)*^g _ ^^^w,B//(Pr)^^ 

(UcVRe)*^g^0 (h’/3);^g^Q 


= <t.(A,B//(Pr)„) 


(D2) 


where c|) is given in equation (15b) taken from reference 1. The second term 
in the numerator may be rearranged as 

(1 “ ^l) 1 / 1 Y _ (1 - Cl) ^^o* Re* / 2 tt 

Cl CRe)„^g - 2 ci Re* (Re)„ g / (ReoXw')* 

w,B=o 

Since Re* = B*/(Pr)^^ and (uc')„ g_Q = /(2 /tt) (ReoX„' ) *, and from equation (C5), 
it follows that 


(1 - Cl) 

ci(Re)w,B 


( ^ ^ 

\ Cl - Cl) , 

f«u'l 

1 j 

fM 

Vup'/Re) 

1 2ci ' 


o 
/ — \ 

o 
' — » 

U j 


(D4) 


The characteristic thicknesses of the present solution (eqs . (B2) and 
(C5)) and the proper boundary conditions require that equation (Dl) should be 
evaluated for the condition 
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since the characteristic length determined from the velocity distribution, for 
u/ug - 0.99, is always the same as the shock-layer thickness. 

A general form of equation (Dl) incorporating equations (D2) through (D5), 
thus, becomes 

I ^ ^ (k/cp)* ; 


In hypersonic flows where e << 1, the factor y is close to lonity. 


thus , 


For T^/Tg - 0, only the first term in equation (D6) need be retained. 




The result, exactly analogous to equation (D6) , can be derived for the 
cases where T^/Tg ^ 0; that is, from equations (B8) and (C7) , equation (D3) 
can be written 


(1 - Cl) / 1 

\Uc'/Re 


(1 - cj) Re** / 2tt 

2ci Re** (Re)^^g / (RCoXw')^ 


w,B=0 



(D8) 


which is the same as equation (D4) . Equation (D6) which is the result shown 
in equations (15), therefore, is applicable for all wall temperature. 


24 


REFERENCES 


1. Yoshikawa, Kenneth K.: Linearized Theory of Stagnation Point Heat and 

Mass Transfer at Hypersonic Speeds. NASA TN D-5246, 1969. 

2. Dorrance, W. H. : Viscous Hypersonic Flow. McGraw-Hill Book Co., 1962, 

3. Howe, John T.; and Sheaffer, Yvonne S.: Mass Addition in the Stagnation 

Region for Velocity up to 50,000 Feet Per Second. NASA TR R-207, 1964. 

4. Anfimov, N. A.; and Aistov, V. V.: Heat Transfer, Friction, and Mass 

Transfer in a Laminar Multicomponent Boundary Layer with Injection of a 
Foreign Gas. International Chemical Engineering, vol. 6, no. 1, Jan. 
1966, pp. 137-144. 

5. Pappas, C. C.; and Lee, G. : Effects of Mass Addition of Various Gases on 

the Heat Transfer and Surface Pressures on a Blunt Cone in Hypersonic 
Flow. AIAA paper 69-716, AIAA Fluid and Plasma Dynamics Conference, 

San Francisco, June 16-18, 1969. 

6. Chen, S. Y.; Baron, J.; and Mobley, R. ; Stagnation Region Gas Injection 

in Low Reynolds Number Hypersonic Flow. Proc. 1967 Heat Transfer and 
Fluid Mechanics Inst., Paul A. Libby, Daniel B. Olfe, and Charles W. 

Van Atta, eds . , Stanford, Calif., Stanford Univ. Press, 1967, pp . 34-57. 

7. Yoshikawa, Kenneth K. : A Simplified Solution of Stagnation Point Heat, 

Shear Stress, and Mass Transfer at Hypersonic Speeds. Paper presented 
at the 6th U. S. National Congress of Applied Mechanics, Harvard Univ., 
Cambridge, Mass., June 15-19, 1970. 


NASA-Langley, 1911 12 A- 36 7 7 


25 



National Aeronautics and Space Administraiion 
Washington, D.C 20546 


official business 


FIRST CLASS MAIL 


04U 001 37 51 3DS 71058 00903 

AIR FORCE WEAPONS LABORATORY /WLOL/ 
KIRTLANO AFB, NEW MEXICO 87117 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE administration 


ATT E. LOU BOWMAN, CHIEF, TECH. LIBRARY 


POSTMASTER: 


If UndeJiverabie (Section 158 
Postal Manual ) Do Not Return 


''The aeronautical and space activities of the United States shall be 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the tvidest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.*^ 

— National Aeronautics and Space Act of 1958 


NASA SOENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Sdemilic aoA 

technical information considered important, 
complete, and a lasting contribution to' existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs^ 
Technology Utilization Reports and 
Technology Surveys, 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



